The Effect of Air on Granular Size Separation in a Vibrated Granular Bed 
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Using high-speed video and magnetic resonance imaging (MRI) we study the motion of a large 
sphere in a vertically vibrated bed of smaller grains. As previously reported we find a non-monotonic 
density dependence of the rise and sink time of the large sphere. We show that air drag causes relative 
motion between the intruder and the bed during the shaking cycle and is ultimately responsible for 
the observed density dependence of the rise time. We investigate in detail how the motion of the 
intruder sphere is influenced by size of the background particles, initial vertical position in the bed, 
ambient pressure and convection. We explain our results in the framework of a simple model and 
find quantitative agreement in key aspects with numerical simulations to the model equations. 

PACS numbers: 45.70.-n, 64.75. +g, 83.80.Fg 
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I. INTRODUCTION 

Granular size separation is one of the most prominent 
phenomena in granular physics 0- 12- 01 ■ When a granular 
mixture is agitated via an external driving force, grains 
of different sizes separate into distinct regions of the con- 
tainer. This behavior sets granular materials apart from 
ordinary fluids which typically mix to increase entropy. 
A key reason is that granular media are far from thermal 
equilibrium, and that therefore the dynamics set up by 
the driving force govern their behavior. The phenomenon 
of granular size separation is well known. It manifests it- 
self in a wide array of granular systems such as chute 
flows, avalanches and rotating drums 0, H H H EJ ■ 

Here we concentrate on a vertically vibrated bed in 
which a large spherical object, the "intruder", is embed- 
ded. The excitation typically causes the intruder to rise 
to the top of the bed, which is commonly referred to 
as the Brazil Nut Effect E3- However, the detailed be- 
havior is more complex. In particular the rate of rising 
is strongly density dependent and in some experiments 
light intruders have been found to sink [lllll3| . Efforts to 
identify the underlying mechanism for the intruder mo- 
tion have been stymied b y th e need to account for this 
wide variety of behaviors |l2j • 

Until recently, the rising or sinking of an intruder has 
been treated as distinct phenomena. Although the in- 
terstitial fluid, typically air, has been shown to strongly 
influence this density dependent behavior of the intruder 
[l3L llH Hil Hr3 | , the details of how it acts and whether 
it is the cause of this behavior have only been re- 
solved recently EH . While some rep orts have suggested 
that air is unimportant El El El El EH, others ac- 
knowledged, its role in determining the intruder motion 
ES El 111 ES El El El, but gave conflicting expla- 
nations. There have also been studies of closely related 
systems wit h multiple i ntruder , quasi 2D and bi-disperse 
systems 
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FIG. 1: Sketch of the setup. H is the total height of the bed 
and h s the initial height of the intruder measured from its 
top. A vacuum gauge measures the background air pressure 
inside the cylinder. 



This paper presents a systematic study of the effect of 
air on size separation in a three-dimensional, single in- 
truder system and expands on a letter published earlier 
EH- I n general, size separation is caused by a variety 
of mechanisms (for a recent review see 0). Our system 
was chosen such that convection is the dominant trans- 
port mechanism in the absence of air. The separation 
mechanism is identified by employing magnetic resonance 
imaging (MRI) and high speed video. This allows us to 
map out in detail the dependence on system parameters 
and to establish a phase diagram that delineates the ris- 
ing and sinking regimes. We then develop a model that 
can explain the key features of our results and also pro- 
vides a unifying framework for seemingly conflicting or 
unconnected observations made in the literature. In par- 
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FIG. 2: Flow fields of the two observed convection patterns: 
(a) axisymmetric, wall-driven convection with a centered heap 
surface. Thin downstream at the walls and an upstream in 
the middle, (b) Asymmetric, air-driven convection with a 
slanted heap surface. Both images are to scale and represent 
flow patterns for shaken beds of d — 0.5 mm glass media. Cell 
radius is R — 4.1 cm and total bed height before shaking is 
H = 8.5 cm. The height ht = 4.5 cm separates the regions of 
up and down flow of the particles at the wall. 



ticular, it shows the sinking and rising to originate from 
the same mechanism. We test this model through a sim- 
ulation. 



II. EXPERIMENTAL METHOD 

All experiments were performed in an acrylic cylinder 
with an inner radius of R = 4.1 cm. It was mounted on 
a VTS 100 electromagnetic shaker driven by a function 
generator that produced sinusoidal excitations ("taps") 
spaced one second apart. Figure ^ shows a sketch of the 
setup. The acceleration was measured in terms of the di- 
mensionless peak to peak acceleration T = a p - p /2g and 
monitored by an accelerometer (PCB, model 353B03) at- 
tached to the vessel. The frequency of the excitations was 
typically fixed at / = 13 Hz although we explored ranges 
between 10 and 30 Hz with no qualitative change in the 
results. We used different types of media to study the 
effect of density: glass (p m — 2.5 g/ml, MoSci Corp.), 
zirconium oxide (p m — 3.8 g/ml, Glen Mills), DVB resin 
(p m = 1.0 g/ml, Supelco), tapioca pearls (p m = 1.2 
g/ml), poppy and raj agar a seeds (p m = 1.0 and 1.2 g/ml, 
respectively). Here, p m is the density of the bed particle 
material The range of sizes we used as our bed medium 
was between 0.25 and 2.0 mm. The intruder particles 
were 25 mm hollow polypropylene spheres (Euro-Matic 
Plastics Inc.) filled with varying amounts of material 
such as lead shot to tune the density. The intruder was 
placed along the cylinder axis at height h s (Fig. 2J. Af- 



ter each run, the cell was emptied and refilled to avoid 
compaction effects. 

Two different cells were used: One had a smooth inner 
wall, the other was roughened by gluing glass beads to 
its surface, thus allowing us to study the effect of wall- 
driven convection on the intruder motion 0, HH • During 
the experiments the cylinder was closed with a lid that 
had a pressure gauge and quick release valve mounted 
on it. The cell could be evacuated to pressures between 
0.13 and 101 kPa. The vacuum pump was then discon- 
nected via the valve to avoid vibrations caused by dan- 
gling tubes. 

A thin plastic straw attached to the intruder that ex- 
tended above the bed surface enabled us to track the 
intruder motion inside the bed. We recorded the motion 
of the straw with a high-speed video camera to measure 
the intruder trajectory during the shaking cycle as well as 
the net displacement of the intruder after each tap. Mea- 
surements of the rise time with and without the straw did 
not show any difference within experimental accuracy.. 

Two different techniques were employed to probe the 
interior of the three-dimensional particle bed. One of 
them involved non-invasive magnetic resonance imaging 
(MRI). The bed was layered with MRI active poppy seeds 
and MRI passive rajagara seeds ^t|- As a result, an MRI 
image of an axial cut through the bed shows a stack of 
bright and dark bands. After each tap and once the sys- 
tem has settled we take a MRI snapshot, imaging a ver- 
tical slice through the center. From the displacement 
of the bands between taps we can then deduce the flow 
field. While the size and density differences between the 
two types of seeds is small, it is not clear a priori that 
they do not matter. The poppy seeds are kidney shaped 
whereas the rajagara ones are more spherical, although 
both are of comparable size (d rs 0.7 mm). In order 
to ensure that our MRI results are unaffected by these 
differences and hold for media other than seeds, we em- 
ployed a second technique to visualize the inside of a bed 
of glass beads. In this technique, initially the bed was 
prepared with equally spaced layers of black and white 
glass beads. The flow field at later times was obtained by 
carefully backfilling the cell with water, freezing it and 
then cutting it along the central axis. 



III. RESULTS AND DISCUSSION 

Convective flow patterns. Since convection plays an 
important role in granular size separation in systems such 
as ours |2^|, we first investigated the convective flow. 
There are two types of convection patterns we observe - 
one is symmetric and the other asymmetric (Fig. |5J) . For 
glass media the symmetric roll is observed with d > 0.5 
mm in the smooth cell and d > 0.35 mm in the rough 
cell. This is the well known wall-driven convection roll 
pill^ that has a thin downstream region near the walls 
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FIG. 3: Normalized rise time T r i sc /T* versus p/pm for dif- 
ferent background densities. T* is the rise time at the peak. 
The size of the media is d = 0.5 mm, h 3 = 4.5 cm, V = 5 and 
/ = 13 Hz, H = 8.5 cm. Symmetric convection: (■), Glass 
(p m = 2.5 g/ml); (a), Zirconium Oxide (p m = 3.8 g/ml) in 
rough cell. Inset: asymmetric convection, (■), DVB resin 
(p,„ = 1.0 g/ml) in smooth cell. 



and a large upstream flow in the center of the bed [Fig. 
IHa)]. Small background media, in which air effects are 
more pronounced due to the lower bed permeability, un- 
dergo the asymmetric convection despite careful leveling 
procedures. We found this asymmetric behavior for all 
glass media with diameters d < 0.5 mm in the smooth 
and d < 0.35 mm in the rough cell. In this case a surface 
instability caused by air flow dominates the wall-driven 
convection and a large convection roll spanning the entire 
diameter of the cell sits on top of a small one [Fig. |2fb)]. 
The asymmetric roll tends to drive bed particles towards 
the wall at which point they either move downward or 
upward depending on their vertical position. 

Density dependent intruder rise time. When an in- 
truder is placed in the bed, it will typically rise upon 
shaking. For sufficiently small and light media in which 
air effects are important, the rate of rising strongly de- 
pends on the relative density between the intruder and 
the bed medium p/p m - 

Despite the two different flow patterns in the bed the 
behavior of the intruder is qualitatively the same in both 
the symmetrical and asymmetrical patterns as can be 
seen in Fig. [3J The plot shows the rise time as a func- 
tion of the intruder density, where both axes have been 
rescaled with the convective rise time T* and the den- 
sity of the bed p mi respectively. The non-monotonic 
curves in Fi g. Bl are ty pical for air-driven size separa- 
tion [ll l lll l[idll7Llllf . Light and heavy intruders rise 
faster than those at intermediate density. 

We define the density at the peak, p* , to be the point 




FIG. 4: Normalized rise time T r i se /T* versus p/p m for differ- 
ent excitation parameters, boundary conditions, interstitial 
gases and media. In all cases convection is symmetric and 
H — 8.5 cm. For clarity, all the curves except (0) have been 
shifted vertically with separation 0.4. The vertical dashed 
lines are at p/p m = 0.41 and 0.63, respectively. Unless stated 
otherwise, we used these standard parameters: Rough Cell, 
R — 4.1 cm, h s — 4.0 cm, d = 0.5 mm glass beads, F = 5 
and / = 13 Hz. (■), standard conditions; (a), T — 3; (o), 
R = 6.0 cm; (♦), poppy seeds d = 0.7 mm; (a), ambient gas 
is Helium (pheii um r "*- / Pair /7); (0), d — 1.0 mm, smooth cell, 
h s — 4.5 cm. 



at which the rise time T r j se is maximized or diverges as 
shown in later plots. In Fig. ^p*/p m ~ 0.5. The curves 
correspond to different background densities of the same 
size (d = 0.5 mm). Both the glass and zirconium oxide 
bed exhibit axisymmetric convection as depicted in Fig. 
Ufa), while the DVB medium does not, despite having 
the same size and thus the same permeability. However, 
due to its lower density it is more susceptible to air effects 
which causes the axisymmetric roll to go unstable and the 
convective flow becomes asymmetric as illustrated in Fig. 
Pl^b). Instead of observing a peak at p* / p m , we find that 
the rise time diverges, since the convective flow drives the 
intruder with p — p* to the wall where it stays for long 
times (sometimes indefinitely). The inset to Fig. [3]shows 
the divergence at approximately the same relative density 
as the peak in the main panel. This means that p* / p m 
does not depend on the symmetry of the convective flow. 

The ratio p* / p m is almost constant over a range of sys- 
tem parameters (Fig. 0J. The density of the interstitial 
gas and the medium, the size and boundary condition of 
the cell and the excitation strength do not significantly 
influence the peak position, p* / p m always lies between 
0.41 and 0.63. This range of peak positions agrees with 
recent results by other groups in similar systems 0, 0] . 
Around the same relative density Yan et al. further- 
more observed a diverging rise time. However, in their 




FIG. 5: Axial cuts through cell for three different density 
intruders. The bed has been cut after it was saturated with 
water and frozen. The bed medium consists of 0.5 mm glass 
beads, interleaved with black-colored layers. Left: sketch of 
original configuration, (a) p/p m = 0.043, 26 taps; (b) p/pm = 
0.5 = p* J Pmi 26 taps; (c) p/p m = 3.3, 6 taps. In (a) and (c) 
the final intruder positions are denoted by a dotted circle. 
Sample (b) was cut with the intruder still inside the bed. 



system the intruder sank below that density. 

Intruder motion and convective flow. Our previous 
MRI experiments established that intruders with density 
p* move with the surrounding bed while light and heavy 
intruders rise faster than convection ^TJ- MRI is, how- 
ever, limited to imaging bed particles that contain traces 
of water or oil, such as seeds. Therefore we checked this 
important assertion more generally for other particles, 
by freezing a water-saturated sample as explained ear- 
lier. Figure [5] shows axial cuts through the bed for three 
different density intruders. The results confirm the MRI 
data for poppy seeds 0. At a density ratio p/p m ~ 0.5, 
the intruder does not move with respect to its surround- 
ings and follows the convective flow. It has not separated 
from the first black layer it was originally sitting on. The 
heavy intruder shows a characteristic wake below it and 
is well above the first black layer. Nearby particles are 
drawn under the sphere. The light intruder is also clearly 
displaced from its original position but does not exhibit 
the pronounced wake shown by the heavy intruder. This 
indicates different mechanisms for the light and heavy in- 
truder by which they rise faster than convection. This be- 
havior also holds for the asymmetrically convecting bed. 
Intruders near pj p m « 0.5 follow the convective flow that 
drives them towards the wall. 

Effects of air and phase diagram. MRI images of in- 
truders rising in an evacuated bed (Fig. |HJ) show that the 
density dependence of the rise time is due to the presence 
of air. Intruders with different densities all rise with the 
surrounding bed. There is no relative motion between the 
intruder and its vicinity. The only transport mechanism 
is global convection [ll HI llH . 

At atmospheric pressure, on the other hand, the air 
enables intruders which are lighter or heavier than p* to 
rise faster than convection (Fig. [SJ). This result originates 
from the interaction between the bed, the intruder and 
the interstitial air which we will discuss later. However, 
this is not the only possible outcome. The same system 
can exhibit sinking behavior. This is one of our central 
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FIG. 6: Absence of density dependence in vacuum. The 
picture are MRI images of axial cuts through the bed for dif- 
ferent density intruders at p = 0.13 kPa. System dimensions: 
R — 3.8 cm, H = 4.2 cm and D — 1.6 cm. The excitation 
parameters were / = 10 Hz and T = 3. Upper row shows ini- 
tial configuration. Lower row shows system after 8 ± 1 taps. 
Column (a) p/p m = 0.05, (b) p/p m = 0.45, (c) p/p m = 2.4. 



results and is summarized in Figs. Hand 03 These data 
show how pressure and initial vertical position in the bed 
determine whether the intruder moves up or down. 

In order to test the influence of convection we per- 
formed the same experiments in rough and smooth cells, 
which differ in convection speed and flow pattern, keep- 
ing all other parameters the same. In both cases we used 
a 0.5 mm glass medium. The rough cell exhibits symmet- 
ric wall-driven convection. In the smooth cell convection 
is significantly reduced and the flow is asymmetric as de- 
picted in figure|2Ib). In Figs. 0and|Hlthe initial starting 
height is just above the middle of the bed at h s = 5.5 
cm. In both cases the curve either peaks or diverges at 
P* I Pm ~ 0.5. As the pressure is lowered p* / ' p m and the 
associated rise times decrease. It is known from Kroll's 
work that, due to air drag, the bed does not lift off 
as high as it would in vacuum. This has a direct con- 
sequence on the wall-driven convective flow. There is 
less shear with the side walls and thus, convection is re- 
duced. As a result intruders rise faster in vacuum than in 
the presence of air. This leads to the shorter rise times 
we observe at low pressures. In both cells the density 
at which the peak occurs decreases significantly at low 
pressures [ri| . 

We also observe a transition between a rising and sink- 
ing regime in the smooth cell. Below a certain pressure 
intruders with densities below p* sink. The resulting 
rise/sink time curve is shown in Fig. IBIb). The data 
resemble those found earlier by Yan et al. When 
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FIG. 7: Rise time curves at different pressures and the phase 
diagram in the rough cell. T — 5, / = 13 Hz and H = 8.5 cm. 
Intruder rise time Hisc versus p/p m at different pressures P 
in the rough cell at h s = 5.5 cm: (■), 101 kPa; (•), 47 kPa; 
(▼), 13 kPa; (a), 6.7 kPa; (0), 0.13 kPa. Inset: Phase dia- 
grams delineating the rising and sinking regimes for d — 0.5 
mm beds at various pressures P. Above each phase bound- 
ary, intruders rise. Shaded area shows the sinking regime at 
ambient pressure. Solid and dashed lines indicate total bed 
height and starting height, respectively. (■), 101 kPa; (o), 27 
kPa and (□), 0.13 KPa. 
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the pressure is lowered even further the intruders rise 
again and the density dependence of T r i sc vanishes. 

This crossover between rising and sinking is not only 
controlled by pressure but also by the vertical position 
of the intruder in the bed. If the intruder is initially 
placed below a certain crossover height h c it sinks. In 
general, h c depends on several other parameters as well: 
The pressure P, the size of the medium d, the density 
p m , the total height H and the wall roughness. 

In order to delineate the rising and sinking regimes, 
we map out phase diagrams in terms of the normalized 
starting height, (h s — D)/D, and the normalized den- 
sity, p/pm (insets to Figs. 0and|Hl and main panel of 
Fig. [5J . The lines connecting the data in these phase di- 
agrams give the experimentally determined, normalized 
crossover height (h c — D)/D. Started above the lines, in- 
truders rise, started at or below the lines they sink. The 
grey areas in the insets of Figs. 0a) and|8£a) denote the 
sinking regimes at atmospheric pressure. As the pressure 
is lowered h c initially increases and then drops. The in- 
truders, regardless of their density, rise from the bottom 
of the cell (h c = D) when P < 0.13 kPa. In this regime 
air effects do not play any role. In the rough cell, h c is 
lower and the pressure dependence not as pronounced as 
in the smooth cell. We believe that strong convection in 
the rough cell gives an upward bias and therefore lower 



FIG. 8: Rise/sink time curves at different pressures and the 
phase diagram in the smooth cell. P = 5, / = 13 Hz and 
H — 8.5 cm. (a) Intruder rise time T r i so versus p/p m at dif- 
ferent pressures P in the smooth cell at h a — 5.5 cm: (□), 101 
kPa; (A), 40 kPa; (o), 27 kPa; (♦), 0.13 kPa. Inset: Phase di- 
agrams delineating the rising and sinking regimes for d = 0.5 
mm beds at various pressures P. Above each phase bound- 
ary, intruders rise. Shaded area shows the sinking regime at 
ambient pressure. Solid and dashed lines indicate total bed 
height and starting height, respectively. (■), 101 kPa; (x), 
27 kPa; (A), 2.7 kPa; (V), 0.67 kPa in the smooth cell, (b) 
The divergent sinking regime. (■), 13.3 kPa with p* — 0.28 
and (□) 2.7 kPa with p* — 0.24. Same conditions as in (a). 
The vertical dashed lines correspond to p* . 



crossover heights compared to the smooth cell. The dif- 
ference between the rough and smooth cell might have 
another cause. Increased shear at the walls could cause 
bed dilation and therefore provide a shortcut for the air 
due to increased permeability. This would lead to a de- 
crease in air effects in the rough cell. 

The sinking of the intruder is caused by pressure gra- 
dients across the bed exerting a net downward pull on 
the intruder 0. Therefore, the phase diagram changes 
with the permeability of the bed, since the magnitude of 
the pressure gradient is inversely proportional to the per- 




FIG. 9: Phase diagram for different sized glass media in the 
smooth cell at P = 101 kPa. T = 5, / = 13 Hz and H = 8.5 
cm. (a), 0.25 mm; (o), 0.5 mm; (■), 1.0 mm; (0), 2.0 mm. 
Inset: Crossover height h c versus total height H normalized 
by D for 0.25 mm glass medium in the smooth cell. 



meability according to Darcy's law. We can adjust the 
permeability k by changing the size of the background 
medium, since k oc d . Figure shows the phase dia- 
gram for glass beads of different sizes. As the medium 
size increases, the sinking regime shrinks. This is ex- 
pected since in the limit of large permeability the pres- 
sure gradients are too weak to induce sinking. Indeed, we 
find that above d = 2 mm all intruders rise with convec- 
tion and air effects become negligible. The inset in Fig. 
El shows (H — h c )/D versus the normalized total height 
H/D. The curve saturates as H increases. 

Size dependence. In the results shown so far, the rela- 
tive size between the intruder and the medium, D/d, has 
been kept at 51. It is vital for understanding the effect 
to see whether the strong density dependence is still ob- 
served as D/d approaches 1. Figure ITUT a) shows the size 
dependence of the rise time for different intruder densi- 
ties at atmospheric pressure. Size affects the rise times 
notably. Both size and density dependencies become less 
pronounced for smaller Did. This might also explain the 
results of Huerta et al. |l8|. They find the convective 
rise time to be faster than most intruder rise times. How- 
ever, they measured convection in the absence of intrud- 
ers which is different from when an intruder is present. 
Fig. HUf a) clearly shows that convection rise time in the 
absence of an intruder (dotted line) is lower than the rise 
time of the intruder at the peak (pi p m = 0.47). From our 
previous MM measurements [ljj and frozen-bed images 
(Fig. 03) we know that this intruder rise time reflects 
the actual convection time. We believe that the large 
sphere blocks air flow, thereby decreasing permeability 
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FIG. 10: Intruder rise time Trisc versus relative size D/d for 
different densities at two different pressures: (a) atmospheric 
pressure P = 101 kPa and (b) low pressure P = 0.13 kPa. 
r = 5, / = 13i?z, H = 8.5 cm and h s = 4.0 cm. The num- 
bers next to the line indicate the relative density. The dotted 
line indicates the convective rise time without the intruder. 
Relative densities in panel b) are: (()), p/p m = 0.043; (V), 
p/pm = 0.26; (A), p/pm = 0.47; (o), p/p m = 0.88. The heav- 
iest intruders (p/ p m = 3.24) stop rising roughly one intruder 
diameter below the surface. 



which in turn slows down convection. At relative sizes 
D/d < 10 the influence of the intruder on the bed dimin- 
ishes and rise time of the intruder at the peak approaches 
the convective rise time without the intruder. Therefore, 
at small D/d air-driven size separation becomes negligi- 
ble and convection dominates [24j. 

At pressures below 0.13 kPa [Fig. HOf b)]. the rise times 
are nearly independent of p/p m and D/d and are close 
to the convective rise time in the absence of an intruder. 
The slight increase in T r i sc at low D/d might be due to ge- 
ometric arching effects as discussed in [27J, |28|, |2|| yfj, |31| . 
Also, the heaviest intruder (p/p m = 3.24) does not sur- 
face at low pressures. It stops rising just below the sur- 
face and appears to sink back into the fluidized layer. At 
atmospheric pressure, the heavy intruders have a longer 
flight time than the bed due to their larger inertia so 
that the bed is already condensed when the sphere hits 
the surface. This explains why it does not sink back in as 
much as it does at low pressure when the upper surface 
is still fluidized at impact. 

Intruder trajectories. In order to elucidate the mecha- 
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FIG. 11: Depth dependence of intruder rise time in 0.5 mm 
glass medium. Total bed height H = 13.5 cm, F = 3.5 and 
/ = 13 Hz. Upper panel (a): T r i sc versus z/D for different 
density intruders (D/d — 50) at atmospheric pressure. (•), 
p/pm = 0.043; (■), p/pm = 0.52 w p*/p m ; (T), p/p m = 3.3. 
Inset shows the velocity at different depths. Lower Panel (b): 
Same plot as above, but at p = 0.13 kPa. (o), p/p m = 0.043; 
(□), p/p m = 0.52 ?s p* / p m \ (V), p/pm = 3.2; (0), convection 
without intruder. Solid line is a linear fit with slope 0.058 
tap - . 



nisms by which light and heavy intruders rise faster than 
convection, we tracked their positions during the ascent. 
Panel HIT a) shows the rise time as a function of depth 
for different density intruders at atmospheric pressure. 
The inset shows the displacement per tap versus depth. 
Deep in the bed the light intruder moves as fast as the in- 
truder at the peak density p* . It then speeds up as it ap- 
proaches the surface. The heavy intruder is always faster, 
but slows down near the surface, because it falls back into 
the fiuidized top layer at the end of the cycle. Below 0.13 
kPa, the curves collapse and the only transport mecha- 
nism is convection [Fig. Illf bl]. This is consistent with 
the convective rise time equation z(t) = £ ln(l + T r j SC /r) 
found by Knight et al. 26] given that in the present ex- 
periments the bed height H is comparable to the depth 
of the convection roll, £. This implies that T r i se <C t, 
leading to z « (£/r)T r i sc . The data in Fig. IllT b) are 




FIG. 12: Intruder displacement in the shaker frame. The 
vertical dashed lines delineate the part I and part II of the 
cycle (see text), (a) Three different densities with h 3 = 7.0 
cm in the rough cell. (□), p/p m = 0.043; (o), p/p m — 0.52 ~ 
P*/Pm', (a), p/p m = 3.3. The solid, dashed and dotted lines 
are simulation results for trajectories of the light, medium and 
heavy intruder (Section V). (b) Light intruder (p/p m = 0.043) 
in smooth cell rising at h a = 6.5 cm > h c , (□), and sinking 
at h s = 3.5 cm < h c , (o). The solid and dashed lines are the 
respective simulation results. The packing fraction was chosen 
to be <f> — 0.59 and <f> = 0.63 for the rough and smooth cell, 
respectively. For the sinking intruder in (b), the simulation 
was stopped at t = 36 ms when the bed volume underneath 
the intruder hits the bottom of the cell. 



well fit with £/r = 1.5 mm/tap (solid line). 

Using high-speed video to image the trajectory of the 
intruder during one shaking cycle we arrive at a more de- 
tailed understanding of the mechanism by which the in- 
truders rise or sink. Figure lT2T a) shows the trajectories of 
intruders for three densities. The vertical lines delineate 
two parts of the shaking cycle. In part I, air flows down 
to the gap that opens up at the bottom the cell [3^, |4l| . 
In the second part, the gap starts to close and air has to 
leave the bed. At p* the intruder has the smallest net 
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FIG. 13: Intruder velocity in the lab frame at: (a) atmo- 
spheric pressure P — 101 kPa and (b) low pressure P = 0.13 
kPa. The vertical dashed lines delineate the part I and part 
II of the period (see text). Three different densities with 
h s — 6.0 cm in the rough cell. (□), p/p m = 0.043; (■), 
p/pm = 0.52 ~ p*/p m ; (A), p/pm = 3.3. The accelerations 
denote the the slopes of the corresponding solid lines. 



upward displacement. The heavy intruder rises higher 
than all the others, as its inertia is large compared to air 
drag. Therefore, its net upwards displacement happens 
in part I. The light intruder on the other hand does not 
rise higher than the one at p* . However, it slows down 
in part II, leading to a higher net upward displacement 
than the intruder at p*. 

The velocity trajectories of the three intruders in the 
lab frame [Fig. I13f a)] show that the heavy intruder expe- 
riences little retardation. The slope of the velocity curve 
is close to — g as denoted by the solid line throughout 
part I and II. It only slows down once the bed starts to 
condense after part II. Conversely, the lighter intruders 
are subject to a significant drag force as evident from 
their high downward acceleration: the first part of the 
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FIG. 14: Trajectory of light intruder (p/p m = 0.043) in the 
lab frame in d — 0.5 mm glass medium at different depths in 
the rough cell. T = 5, f = 13 Hz. (A), h s = 4.0 cm; (□), 
h s = 5.5 cm; (T), h s = 7.5 cm. The black line is a parabola 



with acceleration a 
versus depth. 



-27 m/s . Inset: (■), terminal velocity 



velocity curves is well fit by the — 3g slope of the solid 
black line. They both slow down at the beginning of part 
II when the air flow reverses. 

Figure [r^T b) shows the trajectories for light intruders 
(p <C p*) above and below the crossover height. Both in- 
truders are slowed down in part II by the escaping air. In 
part I, however, the sinking intruder does not rise as high 
as the ascending intruder which causes a net downward 
displacement after part II. This observation is crucial for 
understanding the underlying mechanism. 

Several studies 0, 0, H3, have suggested that 
the right hand side of the peak is caused by inertia. In 
this model heavy intruders penetrate more easily through 
the bed which is held back by shear with the walls. In 
vacuum, where gravity and wall-shear are the only forces 
acting on the bed, our data for the intruder trajectories 
show no evidence of this. All intruders experience the 
same acceleration close to g during flight [Fig. ICT b)]. 
Therefore, the wall shear transmitted through the bed 
exerts only a small force on the intruder in the absence 
of air with little or no density dependence. 

Figure 1141 shows the light intruder trajectory for dif- 
ferent depths in the lab frame. As in Fig. ^1 the light 
intruder slows down in part II of the cycle and reaches 
a terminal velocity. There is a pronounced kink at the 
onset of the slowdown. The closer the intruder is to the 
surface the more it slows down, which increases the up- 
ward displacement and is in agreement with our previous 
observation in Fig. 1111 The inset shows how the terminal 
velocity decreases as the intruder nears the surface. 
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Air drag. In the following we look at the influence of 
the pressure gradients inside the bed on the motion of the 
intruder. These pressure gradients can be significant in 
vibrated beds [39|. To see its effect on the intruder mo- 
tion we performed a simple experiment: Two intruders, 
one with p w p* and the other with p p*, were put on 
top of a bed. After tapping, the trajectories of the two 
intruders differ notably at atmospheric pressure. While 
the heavy intruder is visible throughout the cycle, the 
lighter one is "sucked" into the bed and disappears. At 
the end of the shaking cycle, the light intruder reappears, 
while the heavy one burrows into the bed due to its large 
inertia |43j |. In vacuum, the air drag vanishes and they 
rise to the same height. No relative motion is observed 
between the two intruders (except that at the end of the 
cycle the heavy intruder dives into the bed again). To 
see if Stokes drag plays any appreciable role, the intrud- 
ers were put on a plate and subsequently tapped. They 
rose to the same height within our experimental resolu- 
tion demonstrating that Stokes drag on the intruder is 
negligible. 

Summarizing, we find that the drag on the intruder in- 
side the bed is substantial and originates from the pres- 
sure gradient set up by the bed during the shaking cycle. 
In the absence of the bed, where Stokes drag is applica- 
ble, the air drag has no measurable effect on the intruder. 

Huerta et al. 18] proposed a different model for the 
formation of the peak in the rise time versus pj p m . They 
suggest that the peak results from a competition between 
convection and inertia effects, both of which should also 
occur in vacuum. However, in our experiments the den- 
sity dependence always vanished at low pressures. We 
checked this result explicitly for 2 mm Tapioca pearls, 
a medium that is similar to the 3 mm pearls used by 
Huerta et al. The peak observed at atmospheric pres- 
sure flattened out in vacuum. The reason why air is still 
important for these relatively large and porous media is 
that Tapioca is a very light material (p rn = 1.2 g/ml) and 
much more susceptible to drag forces than glass beads 
(p m — 2.5 g/ml) of comparable size. 

The literature on granular size separation mainly deals 
with two systems: The single intruder on one hand and 
binary mixtures which separate into two phases on the 
other hand. In the former system the interaction is only 
between the intruder and the bed, while in binary mix- 
tures the interactions between the particles that make up 
one species need to be taken into account as well. In a 
model proposed by Hong and coworkers [^3] the separa- 
tion of binary mixtures depends on the granular temper- 
atures of each species. The species with the larger diam- 
eter condenses at the bottom of the container - the so- 
called reverse Brazil nut effect - if D/d ^ Pd/pD, where 
Pd and pa are the densities of the small and large species, 
respectively. This model was shown to be in agreement 
with a recent experiment on separating binary mixtures 
33] with particle diameters in the range 2 — 22 mm. 



In order to test for the onset of interactions between 
multiple intruders in our system we investigated cases 
of up to nine intruders in a variety of initial geometri- 
cal arrangements. All intruders were placed at the same 
starting height in a large, rough cell (R = 6.0 cm) with 
glass beads as bed material. We used d — 0.5 and 1.0 
mm beads and intruder sizes ranging from D / 'd — 12.5 to 
51, and density ratios Pd/pD between 0.3 and 2. There- 
fore the condition derived by Hong et al. for the reverse 
Brazil nut effect to occur was always satisfied. However, 
none of the configurations sank. Moreover, the intrud- 
ers moved as a compound with rise times that did not 
show any appreciable difference to our single intruder ex- 
periment. This held independently of the initial spatial 
arrangement of multiple intruders. This result has two 
implications: If the model proposed by Hong et al. js^] 
were applicable, then many more intruders are needed to 
establish a granular temperature that gives rise to con- 
densation at the bottom. Secondly, in the regime of our 
experiment the size ratio D/d is not the only relevant 
parameter and air effects clearly cannot be neglected. 

IV. MODEL 

The experimentally established connection between 
the density dependent rise times and the presence of pres- 
sure gradients during shaking can be understood using a 
simple model that we develop in this section. This model 
takes into account the interactions between the bed, the 
intruder and the air flow during each shaking cycle and 
enables us to deduce the net displacement of the intruder 
after each tap. In the next section we will incorporate the 
model into a simulation to provide a more quantitative 
comparison to our experimental results. 

Each shaking cycle is divided into two parts: During 
Part I, the bed lifts off the bottom of the cell and a gap 
beneath opens up. The low-pressure region formed at the 
bottom causes air to flow down. As a result the bed and 
the intruder experience a downward force. Smaller me- 
dia with lower permeability increase that air drag. The 
second part of the cycle starts when the gap starts to 
close again, the pressure gradient reverses sign and air 
flows u pwa rds until the bed hits the bottom of the cell 
|39i I40L l4l| . Let us now make some simplifying assump- 
tions: We treat the bed as a porous solid with a constant 
packing fraction </>; Any gap caused by relative motion 
between intruder and bed is immediately filled up; The 
pressure gradients created during the shaking cycle do 
not depend on the position within the bed; Horizontal 
pressure dependence due to the presence of walls is neg- 
ligibly small; Finally, we do not consider convection. 

Within this model the vertical pressure gradient across 
the bed is governed by Darcy's law 
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where k is the permeability of the bed, p the viscosity of 
the interstitial fluid and u the fluid velocity outside the 
bed. The pressure P(z,t) is a function of time and the 
vertical coordinate z only. Using Eq. Q we find the drag 
force on the bed medium F m and intruder Fi nt : 



F m = dP/dz ■ V m = —y———Tnb / ,i 



<ft ■ Pn 



F„„ = (j> PzdS = dP/dz-V int 



dP/dz 



-tnint 



(2) 



(3) 



In the last integral we used the assumption that dP/dz 
is a function of time only. This drag term, Fi nt , has been 
used before to describe the effect of air on a particle in a 
vibrated bed 0, |3f| 01 • Comparing F m and F int , it is 
clear that the bed and the intruder experience the same 
acceleration when pj p m = (ft. 

We now compare this drag force and the viscous drag 
acting on the intruder. The viscous drag in laminar flow 
is given by Stokes formula: 



tokcs 



(4) 



In our setup, r = 0.0125 m, u « A(2irf) = 0.6 m/s 
and mint ranges from 0.5 to 66.7 g, so F stokcs / (m int ■ g) 
ranges from 5.2 • 10~ 4 to 3.9 • 10~ 6 . Therefore, we can 
safely ignore simple Stokes drag on the intruder. The vis- 
cous drag on the bed has been accounted for in Darcy's 
law, where the bed is treated as a continuous porous 
block. The intruder is also subject to the force due to 
pressure gradients as shown in equation ©■ The typi- 
cal gradients in our bed with 0.5 mm glass bead during 
the shaking cycle can be estimated from Darcy's law QJ: 
dP/dz w 25 • 10 3 Pa/m using fj, = 1.8 • 10~ 5 Pa s for air 
and the Carman-Kozeny relation for the permeability of 
a random packed bed of spheres (3(| 
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Therefore Fi nt /(mi n t ■ g) is between 0.31 to 42. Even 
though Stokes drag does not influence the intruder ap- 
preciably, the forces due to pressure gradients are sub- 
stantial. 

It is useful to estimate the effective acceleration acting 
on the bed in part I. The apex of an object's trajectory 
on which a constant drag force Fd acts and an initial 
velocity v is given by h max = v$/(2g ef f), where 
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after substituting (J2J). 

Substituting p m = 2500 kg/m 3 for glass, a typical 
packing fraction (ft = 0.54 and u = 0.6 m/s we obtain 
g e ff — 27 m/s 2 , close to the acceleration of the light 



intruder in Figs. Et a ) an d El until it slows down and 
the kink develops. We expect the intruder at density p* 
to have the same acceleration as well in part I, since it 
moves with the bed. This is indeed the case as shown 
in Fig. I13f a). The intruder has the same acceleration 
before it slows down in part II. The slow down is not as 
abrupt as for the light intruder. 

At density p* / p m = (ft the accelerations of the intruder 
and the bed are the same, so there will be no relative 
motion. In order to illustrate the consequences of air 
drag on different density intruders it is instructive to look 
at the two parts of the cycle for an extremely light and 
heavy intruder. 

Part I: If p 3> p* , the downward acceleration on the 
bed is greater than for the intruder as shown by Nay lor 
et al. [22J. The intruder will push the bed above it. The 
closer the intruder is to the surface the less material it 
has to displace. If p <C p* , the intruder will be pushed 
against the bed below it. Two things can now happen: If 
the intruder is near the bottom of the bed, less material 
has to be displaced and the intruder can sink with respect 
to the bed. If the intruder is located close to the surface, 
the material below will block the downward motion of 
the intruder and it moves with the bed. 

Part II: For intruders with p ^> p* the time of flight 
is longer than for the bed. While the bed collapses the 
intruder is still in the upward motion with respect to the 
bed. The escaping air has little effect on the intruder and 
the void below it is filled up with bed material. By the 
time the trajectory of the intruder has passed its apex the 
bed has almost collapsed and it will have a soft landing 
on the bed resulting in a net upward displacement. The 
light intruder with p <C p* will be accelerated upwards 
with respect to the bed due to its low inertia. If it is near 
the bottom, it cannot displace much of the material above 
it. Above the crossover height, however, it has enough 
inertia to displace the material above and therefore will 
rise with respect to the bed. 

The net displacement of the intruder after one shaking 
cycle is the sum of displacements in part I and II. There- 
fore, intruders with densities different than p* rise faster 
than convection or, for the case of light intruders started 
initially below a crossover height, will sink. In summary 
our model makes the following predictions: 

• At density p* / p m — (ft the intruder does not move 
with respect to the bed. In our simplified model 
without convection it would stay put. In reality, 
however, there is always some convection and the 
intruder will follow whatever background flow there 
is in the system. If the flow is symmetric wall- 
driven convection, then the intruder will surface in 
the middle of the cylinder. In the case of the asym- 
metric convection, the intruder will be driven to the 
side of the container. 

• At all the other densities, the intruder will expe- 
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rience relative motion with respect to the bed due 
to air drag. The intruder will push the bed above 
or below, depending on the relative accelerations 
during the cycle. The dynamics between the in- 
truder and the bed caused by air drag enable light 
and heavy intruders to move faster through the bed 
than those at p* . 

• There exists a crossover height below which light 
intruders with p < p* sink. 

• The density dependence should vanish when the 
relative size D/d becomes small. Even though the 
acceleration due to the drag force on the intruder 
Pint I mint is independent of D, the absolute value 
of its inertia decreases with D 3 . Therefore smaller 
intruders displace less bed than larger ones with 
the same acceleration. Intruders with small size 
ratio D/d rise with convection as is evident in Fig. 
|10fa). This is also the reason why Knight et al. 
|24| did not observe any deviation from convection 
when they measured intruder rise times. All their 
intruders had a relative size D/d < 10. 

• The validity of this model is independent of the 
phase of the sine wave excitation. In general, this 
is true for any excitation in which the acceleration 
exceeds g and therefore causes the bed to lift off. 
Whenever a gap opens up, the air flows downward, 
independent of whether the cell moves up or down 
at that point. Consequently, any excitation can be 
delineated into two phases in which the air flows 
down first and then up. 

A consequence of these considerations is the existence of 
a crossover height. From symmetry arguments one would 
naively expect it to be in the middle of the bed. However, 
convection and fluidization of the upper layers cause an 
upwards bias in the intruder motion. Moreover, the in- 
ertia of the intruder breaks the symmetry. Despite the 
total impulse during the cycle due to drag being close to 
zero, there will still be a net movement. Even if the in- 
truder did not collide with the bed during flight it would 
not return to its initial position. This is due to the differ- 
ent initial velocities at the beginning of part I and part 
II of the cycle. These factors explain why the crossover 
height is usually not in the middle of the bed. 

Despite explaining main features of the separation 
mechanism the model fails to account for some of our 
observations. One is the decrease of p* with decreasing 
pressure (Figs. 0and|HJl. Even though p* / p m < <f>, it is 
unlikely that (j) would decrease to below 0.4. One pos- 
sible explanation is that the pressure gradients are not 
spatially uniform at lower pressures. In this case equa- 
tion Q would no longer hold. 

Another feature of the data the model does not capture 
is the non-monotonic behavior of the crossover height h c 
as a function of pressure (insets of Figs. [7|and[HJ|. 
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FIG. 15: Sketch of the vibrating bed indicating the vari- 
ables used in the model. The dark areas above and below 
the intruder are the cylindrical bed volumes displaced by the 
intruder during the shaking cycle, s is the height of the gap 
between the lower edge of the bed and the cell bottom. 



V. SIMULATION 

In order to test key aspects of our model we performed 
a numerical simulation. In particular the computation of 
the pressure distribution in the bed during one shaking 
cycle sheds light on one main assumption in the model: 
the linear change of pressure with depth at all times. 

The starting point is the Gutman model |38ll39l ] for the 
motion of a vibrated porous bed. It essentially models 
a porous piston in a vibrated cell subject to drag forces 
as the air pulsates through the bed. There is no friction 
between the bed and the walls and the air flow is as- 
sumed to be governed by Darcy's law (JIJ, as before. The 
model also allows for the isothermal compression of the 
gas [38j. Figure El shows a schematic view of the system 
and introduces the relevant variables. 

The pressure P(z, t) obeys the diffusion equation (3(| 
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This equation together with Eq. JIJ are subject to the 
following boundary conditions: At the top of the bed at 
z = 0, the pressure P = Pq, where Pq is the ambient 
pressure. At z = H, the continuity equation holds and 
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where p gas ,Q is the density of the gas at atmospheric 
pressure Pq- Since p gas ~ Pgasfi and s(dp/dt) <C 
Pgas,o(ds/dt) this reduces to 
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Note that this approximation implies that the pressure 
variations are small compared to the absolute pressure. 
In our system, the pressure varies by ~ 4 kPa. Therefore, 
we require P ^> 4 kPa. 

The equation of the bed during the flight is given by 



d 2 (s 
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where x(t) = A sm(2ir ft + 5) describes the oscillatory 
motion of the cell in the lab frame. The bed lifts off 
when the acceleration of the cell reaches —g. Thus, 
5 = arcsin(l/r). We integrate the diffusion equation 
numerically subject to the boundary conditions using the 
Crank-Nicholson scheme • 

The first step in the simulation establishes the pres- 
sure distribution inside the bed during a shaking cycle 
as a function of depth and time. Also the equation of 
motion of the whole bed (|10|l is computed. We neglect 
any pressure redistribution due to the presence of the in- 
truder. This would lead to small corrections that depend 
on the intruder geometry. 

The second step determines the rise and sink time of 
different density intruders at a specified starting height 
h s . The simulation integrates the equation of motion 
of three masses: The bed volumes above and below the 
intruder and the intruder itself (Fig. I15JI . These three 
objects (i — 1,2,3) are subject to drag forces and their 
displacement yi in the shaker frame obeys 
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The simulation computes the net displacement after 
each tap, which is the difference between initial and fi- 
nal position of the intruder plus a small, fixed upwards 
displacement to account for convection. The intruder is 
allowed to displace the upper and lower bed volumes dur- 
ing flight, and after the shaking cycle any gaps are filled 
up to freeze the intruder in place. 

In order for the intruder to sink, there must be a net 
downward displacement from its original position at the 
end of each shaking cycle: The intruder must either pen- 
etrate the material below it or, more likely, some of this 
material must be displaced sideways upon hitting the 
bottom of the cell. In order to avoid introducing more pa- 
rameters and assumptions to the simulation for modeling 
the bed displacement we stop the simulation at this point. 
Nevertheless, we suggest that the region inaccessible to 
the simulation is where sinking is possible according to 
our model. 

Figure ^| shows the results of the simulation for our 
standard system: H = 8.5 cm, F = 5, / = 13 Hz, d — 
0.5 mm, D = 2.5 cm, P — 101 kPa and the convection 
displacement is 0.7 mm/tap. The latter can be obtained 
from the experiments as follows. From Fig. [7] we know 
that at p* the intruder traverses 35 mm in 45 taps. Since 




p / p 

FIG. 16: Simulation results of rise time T r i se versus relative 
density p/ ' p m at two starting heights above and below the 
crossover height: (a) h 3 — 7.0 cm and (b) h 3 = 3.5 cm. 
The parameters are: H = 8.5 cm, / = 13 Hz, V — 5, <j> = 
0.63 (vertical dotted line) and the convection displacement 
is 0.7 cm / tap. The inset in (a) shows the phase diagram 
of this system at atmospheric pressure. Solid and dashed 
lines indicate total bed height and the two starting heights in 
(a) and (b), respectively. The gray area indicates the region 
inaccessible to the simulation (see text). 



the convection speed does not change appreciably with 
depth (Fig. 1111) . the displacement due to convection is 
~ 35/45 = 0.7 mm/tap. The packing fraction was chosen 
to be close to the random closed pack limit, (j> = 0.63, but 
any value between 0.5 — 0.63 yields similar results. 

Comparing the simulation (Fig. 116(1 to the experimen- 
tal data from Figs. 0|HI yields quantitative agreement of 
key features: Both the peak and the divergence in Fig. ED 
occur at p* / p m = 0.5 and 0.6, respectively. Even though 
sink times cannot be obtained in a quantitative fashion 
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from the simulation without introducing additional pa- 
rameters, the grey area in Fig. ll6T b^ is consistent with 
the sinking regime from the experiment. 

Moreover, the simulated intruder trajectories closely 
match the high-speed measurements in Fig. I12f a): 
Denser intruders rise higher and the apex occurs at later 
times. Light intruders on the other hand, slow down in 
part II of the cycle. This slow-down is responsible for 
the net displacement of the light intruder with respect 
to the surrounding bed. The trajectories of the light in- 
truders cross those at p* at approximately t = 0.05 s, 
which agrees well with our simulation result. The kink 
in the simulated light intruder trajectory is caused by 
the intruder traversing the gap created above it in part 
I. In the experiment, this gap is likely to get partially 
filled thereby impeding the intruder motion. The simu- 
lated trajectories of the light intruder above and below 
the crossover height [Fig. 112( b)] agree well with the data 
apart from the kink due to the presence of the gap. The 
trajectory of the sinking intruder is only simulated until 
the lower bed volume hits the bottom of the cell. 

The simulation of a shaking cycle stops when the en- 
tire bed hits the bottom of the cell. In the experiments, 
however, the intruder keeps moving past that point as 
shown in Fig. ^] This is because the bed dilates dur- 
ing flight and condenses as the bottom layers hit the cell 
bottom. This condensation takes approximately 0.02 s. 
The heavy intruder keeps moving even longer, since it is 
penetrating the bed due to its large inertia. These ef- 
fects are not considered in our simulation and, given our 
results, they seem to be second order effects. 

The top boundary of the grey sinking regime is not flat 
[inset Fig. HBT a)]. Instead it dips down at low densities. 
Experimentally, however, we do not observe this for the 
0.5 mm glass medium. Nevertheless, Fig. Elsuggests that 
this dip may exist for large permeabilities. 

Due to the restrictions on Pq by virtue of the approx- 
imations made in equation 0, we cannot probe very 
low pressures in the simulation. Down to P = 40 kPa 
P* I Pm only shifted to 0.45, which is a smaller reduction 
than observed in the experiment |l7||. Even though the 
pressure change does not seem to affect p* significantly, 
the crossover height goes up at lower pressures. More- 
over, making convection stronger by choosing a higher 
value for the convection displacement in the simulation 
decreased the crossover height. 

Summarizing our simulation results, we find that im- 
portant parameters, such as peak and divergence posi- 
tion agree well with experiments. The reason why the 
peak does not occur exactly at 4> is due to the fact that 
the pressure is not changing linearly with depth at all 
times. This is a key assumption in our model. The sim- 
ulation shows departures from this assumption. The in- 
truder trajectories are in good quantitative agreement 
with the experimental data (Fig. I12|) within the limits 
of our simulation. We do not expect a full quantita- 



tive correspondence with the experiment due to several 
factors the simulation lacks. It does not account for pen- 
etration of the intruder through the bed, it lacks any 
fluidization during flight and the experimental excitation 
is not perfectly sinusoidal as assumed in the simulation. 
Moreover, the bed masses that are displaced are assumed 
to be columns, which is another idealization. Neverthe- 
less, given the simplifying assumptions, our model cap- 
tures the essential features of the Brazil Nut phenomenon 
and qualitatively reproduces the experimentally observed 
phase diagrams in Figs. [7]and[S] 

A non-monotonic rise time versus density curve that 
peaks at p*/p m ~ 4> has previously also been seen in 
a model proposed by Rhodes et al [lq . These authors 
use an equation of motion for the intruder that includes 
three additional forces: two types of bed retardation 
and an added mass term. It requires three fitting pa- 
rameters that are not directly obtainable from exper- 
iment. Our model on the other hand needs only one 
directly-measurable parameter, the convection displace- 
ment, which has a clear physical origin. 

VI. CONCLUSIONS 

The role of air in vibration-induced granular size sep- 
aration can be dramatic. In media which are sufficiently 
small and light, air effects dominate the process of size 
separation and can cause intruders to sink that would 
rise otherwise. This process is essentially decoupled from 
background convection. MRI, freezing the bed with wa- 
ter, and high speed video measurements enabled us to 
study the flows inside the bed and elucidate the separa- 
tion mechanism. 

A key finding is that even though Stokes drag on the 
intruder in the absence of the bed is negligible, the drag 
caused by bed- induced pressure gradients is not. For one- 
inch diameter intruders in half-millimeter bed material 
the gradients can lead to forces several tens of times the 
intruders own weight. These forces cause relative motion 
with respect to the bed during one shaking cycle. This 
leads to a phase diagram that separates rising and sink- 
ing behavior as a function of depth and pressure. Intrud- 
ers with p/ p m < 4> si n k if started below some crossover 
height. Our model gives a unified description of the ris- 
ing and sinking of an intruder that agrees with previous 
experimental results. 

The air drag also causes a size and density dependence 
of the speed with which intruders rise or sink. We find 
that at p* I p m ~ <f> the intruder moves with the convec- 
tive background. Heavier intruders rise faster than con- 
vection, while lighter intruders rise or sink with respect to 
convection depending on their initial vertical position in 
the bed. This result as well as measurements of the flight 
trajectories of the intruder during a tap are in quantita- 
tive agreement with the model we developed. 
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In this model, heavy intruders rise because bed parti- 
cles are held back more strongly by drag due to pressure 
gradients. Light intruders, on the other hand, are more 
easily buffeted by the air currents flowing through the 
bed; if placed initially near the bottom, this can push 
the intruder downward, whereas if placed near the top 
the intruder can get driven toward the upper surface. 
However, we cannot explain the decrease of p* as the 
pressure is lowered in the framework of our model. This 
might be due to simplifying assumptions, such as con- 
stant packing fraction, lack of fluidization and intruder 
penetration. Despite the fact that our model cannot ac- 
count for all the rich details this experiment yields, it 
does unify previously disjointed pieces of data, mainly 
the rising and sinking of light intruders. 

The air driven size separation is a non-equilibrium phe- 
nomenon. The different dynamics during the two parts 
of the shaking cycle explain the observed results. Even 
though our studies concentrated on the single intruder 
situation, similar considerations can be made for multi 
intruders and binary mixtures where the size difference 
is large and the bed permeability is small. 
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